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Glueball masses in U(l) LGT using the multi-level algorithm 
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a Max-Planck-Institut fiir Physik, Fdhringer Ring 6, D-80805, Miinchen 

The multi-level algorithm allows, at least for pure gauge theories, reliable measurement of exponentially small 
expectation values. The implementation of the algorithm depends strongly on the observable one wants to 
measure. Here we report measurement of glueball masses using the multi-level algorithm in 4 dimensional compact 
U(l) theory as a case study. 
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1. Introduction 

Compact U(l) lattice gauge theory exists in 
two phases. A confining phase at strong coupling 
and a deconfined phase at weak coupling. Monte 
Carlo simulations have established that for the 
Wilson action, the phase transition point corre- 
sponds to (3 = 1.011128(11) [I]. The order of 
the transition has been debated for a long time. 
Recent investigations including finite size scaling 
analysis suggest a weak first order transition . 

While accurate measurements of the glueball 
mass can throw light on the order of the transi- 
tion, such measurements are difficult to perform 
using conventional methods. Glueball masses in 
4d compact U(l) lattice gauge theory were first 
measured by Berg and Panagiotakopoulos [2] us- 
ing correlations between Wilson loops. However 
they could only go to a separation of 2 in lattice 
units for the correlators. In fact except for val- 
ues of the coupling close to the phase transition 
point 3 (where the glueball is lighter), measure- 
ments have been carried out only for small tem- 
poral lengths of the correlators. 

Recently Liischer and Weisz have proposed an 
exponential noise reduction method which ex- 
ploits the local nature of the action and existence 
of a positive definite transfer matrix 0]. In this 
work we apply this method, which lets us go to 
large temporal separations, to glueball correlators 
and obtain results which have very little contam- 
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ination from higher states. A similar study was 
also carried out in [5]. 

2. Glueball correlators in the multi-level 

scheme 

Glueball correlators generically look like 

(C(t,t )) conn = (0(f)O(to)> - (0(t))<0(t o )>. (1) 

The zero momentum scalar and axial- vector cor- 
relators are given by the operators 0\ and O2 
such that 

Oi(t) 

o 2 (t) 



E Re(P y (f,t)) 

x ij=l,2,3 

Vlm(P 4J (x,t)), 



(2) 
(3) 



where Pij is the plaquette in the ij plane. Glue- 
balls with definite momenta k arc created by 



0(M) 



o(x, ty 



(4) 



As a function of the time separation At = t — to 
the glueball correlator is expected to behave like 



<C(t,t„)>c 



iAt 



_|_ e -m(N t -At) 



(5) 



where N t is the extent of the lattice in the time 
direction. Fitting the measured correlator to this 
form, one can obtain the effective mass m of the 
glueball. For zero momentum, the effective mass 
is equal to the rest mass while for non-zero mo- 
mentum, one has to take into account the mo- 
mentum contribution (mfe) to the effective mass 
to obtain the rest mass. 
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In the multi-level scheme, (0(t)0(t n )) is esti- 
mated by ([0(f)] [O(^o)]) where [■ ■ •] denotes an 
intermediate level of averaging called the sub- 
lattice average 4. This scheme requires partial 
updates of the lattice. In contrast to a full up- 
date where all links are updated, a partial update 
affects links only in a part of the lattice with the 
boundary of this part held fixed. The sub-lattice 
averaging reduces the fluctuation of each individ- 
ual operator O to a great extent. This is efficient 
because the small expectation values are now gen- 
erated by multiplication rather than fine cancel- 
lation of positive and negative values of the same 
order. 

As long as (O) = 0, which is true for the axial- 
vector correlator, this procedure works quite well. 
However in the scalar channel, where (O) ^ 
much of this advantage is lost as each expectation 
value is a number 0(1), but the connected part is 
several orders of magnitude smaller than the full 
correlator. To get around this problem, one can 
take the derivative of the correlator in the scalar 
channel to remove the VEV of the plaquettes. So 
let us now take the derivative of the correlator 
at both t and to- 4 Taking dt to be the forward 
derivative and <9 t * o to be the backward derivative 
on the lattice, we get, 



forward derivative 



d t dl(C(t,t )) « -Q 



,-m(t-t ) 



(1 



+e 



-m(N t -(t-t )) r e m 
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Now we are far better suited to apply the multi- 
level scheme. We can now measure 

d t d; o (c(t,t )) = (J2 [Pij(s,t+i) - Pij(x,t)) 

x [Pij&h) - Pij{2, to -1)]). (7) 

S,i,j 

The derivatives were estimated in the sub- 
lattice updates by taking the difference of the 
value of the operator on the updated slice with 
the value of the operator on the boundary. As 
shown in Figure 1, to get the forward derivative 
at t, we used the fixed boundary at (t + 1) and for 

4 In principle one derivative is enough, but in practice we 
found that the efficiency of the algorithm is higher for the 
double derivative compared to the single one. 
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Figure 1. Evaluation of the derivative of the 
glueball correlator. The thick lines are time slices 
held fixed during the sub-lattice averaging. 



the backward derivative, the boundary at (t — 1). 
To get the correlators one has to use two such 
slices (e.g. t and to in Figure 1). In practice we 
held every alternate layer of spatial links fixed 
and estimated the correlators for various tempo- 
ral separations in the same sweep. The only draw- 
back at the moment seems to be the fact that we 
have to consider a minimum separation of two in 
the temporal direction. 

The number of sub-lattice updates is an op- 
timization parameter of the algorithm that has 
to be tuned for efficient performance. This is a 
function of (3. In the range of (3 we looked at, 
we found that 10 to 50 sub-lattice updates were 
sufficient. To compare this procedure with the 
naive algorithm we measured the percentage er- 
ror on the correlators at a value of (t — to) where 
both methods gave non-zero signals. In a similar 
amount of computer time, the multi-level algo- 
rithm produced errors which were about two or- 
ders of magnitude lower than the naive method. 
In Figure [2] we show the %error for a given CPU 
time as a function of sub-lattice updates for the 
scalar and axial- vector channels. 

3. Results 

In Tables ^ and we present our results for the 
zero momentum scalar and axial-vector glueball 
masses. 

The (k ^ 0) axial-vector correlator is sensitive 
to a correlation between two 1-forms. In the de- 
confining region we expect this part to yield in- 
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Figure 2. %error on the scalar (top) and axial- 
vector (bottom) correlator at At — 2 after a lOhr. 
run on a 1.5GHz AMD Athlon PC. 12 4 lattice. 
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Figure 3. Consolidated plot of glueball masses. 
▲ and T : scalar channel (k — 0) from dd*(C) 
and d*d(C). • : axial-vector channel (k — 0). □ 
: photon from the axial- vector correlator. 
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formation about the photon. Indeed the effective 
masses from this correlator are very small in the 
deconfined phase and do not show any variation 
with p. To obtain the rest mass, we computed 
m 2 — m\ assuming the free field dispersion rela- 
tion m 2 , = Y^i (2 — 2cosfci). Within statistical 
errors the rest mass turns out to be zero and this 
we believe is strong evidence for the photon. For 
a more elaborate discussion of our results see 0. 

In Figure [21 we present all the masses along 
with the region where the phase transition is ex- 
pected to take place. 
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Table 1 



Zero momentum scalar glueball masses 



& 


d;d ta (c(t,t )) 


d t dt n (C(t,t )) 


0.990 


1.195 (55) 


1.085 (50) 


1.000 


0.875 (41) 


0.812 (37) 


1.005 


0.682 (29) 


0.693 (29) 


1.010 


0.405 (22) 


0.410 (21) 



Table 2 

Zero-momentum axial- vector glueball masses 
Confining regime 





1.000 


1.005 


1.010 


ma 


1.31 (1) 


1.096 (9) 


0.757 (20) 


Deconfining regime 


a 


1.012 


1.015 


1.020 


ma 


1.068 (38) 


1.279 (24) 


1.47 (8) 



